Jet broadening in unstable non-Abelian plasmas 
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We perform numerical simulations of the SU(2) Boltzmann-Vlasov equation including both hard 
elastic particle collisions and soft interactions mediated by classical Yang-Mills fields. Using this 
technique we calculate the momentum-space broadening of high-energy jets in real-time for both 
locally isotropic and anisotropic plasmas. In both cases we introduce a separation scale which sepa- 
rates hard and soft interactions and demonstrate that our results for jet broadening are independent 
of the precise separation scale chosen. For an isotropic plasma this allows us to calculate the jet 
transport coefficient q including hard and soft non-equilibrium dynamics. For an anisotropic plasma 
the jet transport coefficient becomes a tensor with §l q±- We find that for weakly-coupled 
anisotropic plasmas the fields develop unstable modes, forming configurations where B± > E± and 
E z > B z which lead to > q±. We study whether the effect is strong enough to explain the 
experimental observation that high-energy jets traversing the plasma perpendicular to the beam 
axis experience much stronger broadening in rapidity, Ar/, than in azimuth, A<j>. 

PACS numbers: 12.38.-t, 12.38.Mh, 24.85.+p 



INTRODUCTION 

High transverse momentum jets produced in heavy-ion 
collisions represent a valuable tool for studies of the prop- 
erties of the hot parton plasma produced in the central 
rapidity region [l[ . This is due to the fact that jets couple 
to the plasma causing the jet to broaden in momentum 
space and to lose energy. The magnitude of momentum- 
space broadening and energy loss experienced by a parton 
depends on whether or not one assumes the matter to be 
hadronic or partonic in nature. Hence, it is one of the pri- 
mary observables to ascertain experimentally whether or 
not the plasma has been produced. At very high energies 
it is expected that hard bremsstrahlung processes dom- 
inate the light quark or gluon energy loss however, 
at intermediate energies the inclusion of both collisional 
and radiative processes is necessary in order to make phe- 
nomenological predictions. Here we present first results 
from real-time solution of the SU(2) Boltzmann-Vlasov 
equation for locally isotropic and anisotropic plasmas 
which include both hard (scattering) and soft (classical 
field) processes. 

We first demonstrate that in a locally isotropic plasma, 
one can obtain a cutoff-independent transport coefficient 
q, which measures the square of the transverse momen- 
tum transfer per mean free path. This measurement lays 
the groundwork for determining other jet transport prop- 
erties like the energy-loss spectrum from real-time simu- 
lations. We then perform a similar measurement in a lo- 
cally anisotropic plasma and demonstrate that for a non- 
Abelian plasma unstable modes can cause asymmetric 
broadening of jets. This may be relevant for recent mea- 
surements of di-hadron correlations which provide evi- 



dence for an asymmetric broadening of jet profiles in the 
plane of pseudorapidity (77) and azimuthal angle (0), with 
A77 > A<fi, which has been called "the ridge" by experi- 
mentalists [l|, 0] . Here we will show that this asymmetry 
could partly be caused by unstable plasma modes which 
are induced by the longitudinal expansion of the plasma. 

At the earliest times after an ultrarelativistic heavy- 
ion collision, before thermalization and hydrodynamic 
expansion, the plasma undergoes rapid longitudinal ex- 
pansion. This expansion can lead to an oblate anisotropic 
(ip z ) <C (Pj_)) momentum distribution in the local rest 
frame Ij, [5j. It has been shown that instabilities de- 
velop [5j in such anisotropic plasmas which lead to 
the formation of long-wavelength chromo-magnetic and 
chromo-electric fields. These fields can then affect the 
propagation of hard jets and their induced hard radiation 
field. In the Abelian case soft transverse magnetic fields, 
B±_ , dominate all other field components and therefore in 
this case one expects the longitudinal pressure induced 
by unstables modes to be larger than the transverse pres- 
sure @, 0] ■ This pressure asymmetry causes asymmetric 
broadening of jets with larger broadening along the lon- 
gitudinal or r\ direction. 

The situation is more complicated in non-Abelian plas- 
mas since then, in addition to generating large coherent 
B± and E± domains, one also generates large-amplitude 
longitudinal fields B z and E z . It is therefore not obvious 
a priori that the pressure generated by other field com- 
ponents will result in an asymmetric broadening of the 
jet. Through our numerical simulations we find that for 
oblate parton momentum distributions, that at different 
times cither E z > B z or B± > E± with the net effect be- 
ing a factor of 1.5 stronger longitudinal than transverse 
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broadening. 



BOLTZMANN-VLASOV EQUATION FOR 
NON-ABELIAN GAUGE THEORIES 



We solve the classical transport equation for hard glu- 
ons with SU(2) color charge q — q a t a Q, with the color 
generators t a , including hard binary collisions 

% + gg a F^d; + 9 r bc Al(x)q c d q a] f = C, (1) 

where / = f(x,p,q) denotes the single-particle phase 
space distribution. It is coupled self-consistently to the 
Yang-Mills equation for the soft gluon fields, 



V" = f = 9 J dq qv» f(x,p,q), 



(2) 



with v 1 - 1 — (l,p/p). When the phase-space density is 
parametrically small, / = 0(1), the collision term is 
given by 
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with J = J (2n)32E- • ^ ne ma trix element M. includes 
all gg — > gg tree-level diagrams and color factors as ap- 
propriate for the SU(2) gauge group. 

We employ the test particle method to replace the con- 
tinuous distribution f(x,p,q) by a large number of test 
particles, which leads to Wong's equations Q 

*i(i)=Vi(t), (4) 
Pi(i) = 9q?(t) (E a (i) + Vi (t) x B a (t)) , (5) 
q i (t) = -igv?(t)[A lt (t),q i (t)], ( 6 ) 

for the z-th test particle, whose coordinates are Xj(£), 
Pi(t), and qf(t). The time evolution of the Yang-Mills 
field is determined by the standard Hamiltonian method 
[13 in A = gauge. Our lattices have periodic bound- 
ary conditions and a lattice spacing a, whose physical 
value is fixed by the length of the lattice L — aN s . Di- 
mensionless lattice variables scale such that when the 
number N s of lattice sites is changed, L remains fixed; 
see [HldH f° r details. 

The theory without collisions as given by equations ^ 
coupled to the lattice Yang-Mills equations was first 
solved in ref. (l3j to study Chern-Simons number diffu- 
sion in non-Abclian gauge theories at finite temperature. 
It was applied later also to the problem of g auge-field 
instabilities in anisotropic SU(2) plasmas [111. Il2j. Our 
numerical implementation is based on the improved for- 
mulation detailed in ref. [l2| where the non-Abelian cur- 
rents generated by the hard particle modes on the lattice 



sites are "smeared" in time. This technique makes sim- 
ulations in three dimensions on large lattices possible in 
practice. 

In this paper we go beyond those earlier simulations by 
accounting also for hard (short-distance) collisions among 
particles. The collision term is incorporated using the 
stochastic method [13]. Scattering processes are deter- 
mined by sampling possible transitions according to the 
collision rate in a lattice cell: 
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(7) 



with v rc \ = s/(2E\E2), where s is the invariant mass of 
a gluon pair. The total cross section is given by 



02^2 = 



~T~2~dq 

k ,2 dq z 



(8) 



The momentum transfer is determined in the center of 
mass frame of the two colliding particles from the prob- 
ability distribution 



V(q 2 



1 da 

a 2 ^2 dq 2 



(9) 



In Eq. ([8} we have introduced an infrared cutoff k* for 
point-like binary collisions. To avoid double-counting, 
this cutoff should be on the order of the hardest field 
mode that can be represented on the given lattice, k* ~ 
■k / a. Momentum transfers below k* are mediated by the 
classical Yang-Mills field; a soft scattering corresponds to 
deflection of a particle in the field of the other (s). Note, 
that we use the color averaged expression for the collision 
term. The color charge of a particle is hence not affected 
by a hard collision. 

Physically, the separation scale k* should be suffi- 
ciently small so that the soft field modes below k* are 
highly occupied (lQ |. On the other hand, k* should be 
sufficiently large to ensure that hard modes can be repre- 
sented by particles and that collisions are described cor- 
rectly by Eq. ([3]). Furthermore, unstable modes arise in 
anisotropic plasmas (see below), all of which should be 
located below k*. Since g ~ 1, in practice, we choose 
k* = \/3ir I 'a to be on the order of the temperature for 
isotropic systems, and on the order of the hard transverse 
momentum scale for anisotropic plasmas. Independently, 
one should have m^L 3> 1 and moo a <C 1: the first con- 
dition ensures that the relevant soft modes actually fit on 
the lattice while the latter corresponds to the continuum 
limit. Here, m,^ denotes the soft scale and is given by 



rr4 = g 2 N c 



d 3 P /(p) 
(2tt) 3 Ipl 



g 2 N c — 
Ph 



(10) 



where N c — 2 is the number of colors and n g denotes 
the number density of hard gluons, summed over two 
helicities and N 2 — 1 colors. Also, ph w 3T is the typical 
momentum of a hard particle from the medium. 
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As we have argued above, we shall choose the inverse 
lattice spacing to be on the order of the temperature of 
the medium. Thus, with (jTUJ) the continuum condition 
moo a ^ 1 roughly translates into 



l 2 N r 



T 3 



< 1 



(11) 



In order to satisfy this relation with g ~ 1, in our nu- 
merical simulations below we shall assume an extremely 
hot medium, T 3 » n g . However, this should be viewed 
simply as a numerical procedure which ensures that the 
simulations are carried out near the continuum (or weak- 
coupling) limit. We shall verify below that transverse mo- 
mentum broadening of a high-energy jet passing through 
a thermal medium is independent of T if the density 
and the ratio of jet momentum to temperature is fixed; 
compare to Eq. (fT3")) below. One may therefore obtain 
a useful "weak-coupling" estimate of (resp. for the 
related transport coefficient q) by extrapolating our mea- 
surements down to realistic temperatures. 



JET BROADENING IN AN ISOTROPIC PLASMA 



We first consider a heat-bath of particles with a den- 
sity of n g = 10/fm 3 and an average particle momentum 
of 3T —12 GeV. The rather extreme "temperature" is 
chosen to satisfy the above conditions on N s = 32 • • • 128 
lattices, assuming L = 15 fm. For a given lattice (resp. 
k*) we take the initial energy density of the thermalized 
fields to be 



d 3 k 



kf Boso (k)e{k* -k) 



(12) 



where / B osc(fc) = n g /{2T 3 C{Z))/{e k/T - 1) is a Bose dis- 
tribution normalized to the assumed particle density n g , 
and £ is the Riemann zeta function. This is equivalent 
to the energy density of Bose-distributed particles with 
momenta below the separation momentum k*. The ini- 
tial spectrum is fixed to Coulomb gauge and Ai ~ 1/fc (in 
continuum notation); also, for simplicity we set Ei = at 
the initial time but electric fields build up quickly within 
just a few time steps. 

We then measure the momentum broadening (p\)(t) 
of high-energy test particles (p/3T ss 5) passing through 
this medium. Fig. Q] shows that in the collisionless case, 
C = 0, the broadening is stronger on larger lattices, 
which accommodate harder field modes. However, Fig. [2] 
demonstrates that collisions with momentum exchange 
larger than k* (a) compensate for this growth and lead to 
approximately lattice-spacing independent results even 
when k* varies by a factor of four. 

Figs. Q] and [2] show that the relative contributions to 
(p^_) from soft and hard exchanges can depend signifi- 
cantly on k* , even forp/fc* = O(10). It is clear, therefore, 
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FIG. 1: (Color online) Momentum diffusion caused by 
particle-field interactions only. Additional high-momentum 
modes on larger lattices cause stronger momentum broaden- 
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FIG. 2: (Color online) Momentum diffusion by both particle- 
field and direct particle-particle interactions. The result is 
independent of the separation scale k* . Same parameters as 
in Fig. rrj 



that transport coefficients obtained in the leading loga- 
rithmic (LL) approximation from the pure Boltzmann 
approach (without soft fields) will be rather sensitive to 
the infrared cutoff k* . Fitting the difference of Fig. O and 
Fig. [T] (i.e., the hard contribution) to the LL formula 



d(p] 



df 



C f 8tt' 



log C 



k*' 



(13) 



gives C ~ 0.43, 0.41, 0.31 for k*/T = 2^3, 0.5^3, 
respectively. For the full calculation C ~ 0.61 k*/(s/3T). 

A related and frequently used transport coefficient is 
q It is the typical momentum transfer (squared) per 
collision divided by the mean-free path, which is nothing 
but {p\){t)/t. From Fig.|2J q ~ 2.2 GeV 2 /fm for iV c = 2, 
n g = 10/fm 3 and p/(3T) » 5. Our cut-off independent 
value for q is in the range extracted from phenomcnolog- 
ical analyses of jet-quenching data from RHIC (is] ]. 

In FigJ3]we show that q is indeed largely independent of 
the separation scale k*. In these simulations, test parti- 
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FIG. 3: <J as a function of k* at fixed n g = 5/fm 3 and 
p/(3T) = 16. 
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FIG. 4: q as a function of T at fixed n g 
P/(3T) = 5. 
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JET BROADENING IN AN UNSTABLE PLASMA 

In heavy-ion collisions, locally anisotropic momentum 
distributions can emerge due to the longitudinal expan- 
sion. Such anisotropics generically give rise to insta- 
bilities see 11, 12| for simulations of unstable non- 
Abelian plasmas within the present "Wong- Yang-Mills" 
approach. Here, we investigate their effect on the mo- 
mentum broadening of jets, including the effect of colli- 
sions. The initial momentum distribution for the hard 
plasma gluons is taken to be 

/(P) = %( — J S(p z )exp(-px/ph) , (14) 



with p± = yPx + Py ■ This represents a quasi-thermal 
distribution in two dimensions with average momentum 
= 2ph- We initialize small- amplitude fields sampled from 
a Gaussian distribution and set k* « Ph, for the reasons 
alluded to above. The band of unstable modes is located 
below k*. 

We find that binary collisions among hard particles 
reduce the growth rate of unstable field modes, in agree- 
ment with expectations [161 ] . However, for ph = 16 GeV, 
L = 5 fm, n g — 10/fm 3 , g = 2, TOqo ~ 1/fm and k* 
1.7ph, the reduction of the growth rate is only approx- 
imately 5%, increasing to about 15% when k* w 0.9ph- 
This is due to fewer available field modes and more ran- 
domizing collisions. 

Next, we add additional high momentum particles with 
p x = 12 ph and p x = 6ph, respectively, to investigate the 
broadening in the y and z directions via the variances 
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(15) 



cles were explicitly bunched into colorless jets, such that 
radiative energy loss does not contribute. This explains 
why one does not need bremsstrahlung processes in the 
collision term to obtain a cutoff independent result. Note 
that the magnitude of q ~ 1.3 GeV 2 /fm at n g = 5/fm 3 
and p/(3T) = 16 is smaller than the extrapolation ob- 
tained from Eq. (TlB"]) with the constant under the loga- 
rithm fixed from the previous runs at n g = 10/fm 3 and 
p/(3T) = 5; hence, C effectively depends on the density 
and the jet momentum. 

In Fig. [H we verify that q does not depend on the tem- 
perature T so long as the particle density n g and the 
ratio of jet momentum to temperature, p/T, is fixed. 
Thus, our measurement of q ~ 2.2 GeV 2 /fm may be 
considered as a weak-coupling extrapolation to realistic 
temperatures of T w 300 MeV and jet momenta of about 
4.5 GeV. 



The quantity y 'qh/qi. can be roughly associated with the 
ratio of jet correlation widths in azimuth and rapidity: 
\fIiJV. ~ (Ar;)/(Ao). 

Fig. [5] shows the time evolution of (p\) and of (p 2 z ). 
The strong growth of the soft fields sets in at about 
t ~ 10 rn" 1 and saturates around t ~ 25 mZs due to the 



finite lattice spacing (also see [12j). Outside the above 
time interval the ratio <?l/c?± ~ 1- During the period of 
instability, however, 



3-L 



(16) 



for both jet energies shown in Fig. \5\ We find approxi- 
mately the same ratio for denser plasmas (n g = 20/fm 3 
and n g = 40/fm 3 ). Reducing the number of lattice sites 
and scaling p^ down to 8 GeV gives c?l/q_l ~ 2.1. How- 
ever, these latter runs are rather far from the continuum 
limit and lattice artifacts are significant [T3 ] . 

The explanation for the larger broadening along the 
beam axis is as follows. In the Abelian case the insta- 
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FIG. 5: (Color online) Momentum broadening of a jet in 
the directions transverse to its initial momentum. p z is di- 
rected along the beam axis, p± is transverse to the beam. 
Anisotropic plasma, 64 3 lattice. 



bility generates predominantly transverse magnetic fields 
which deflect the particles in the z-direction [7( . 
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FIG. 6: (Color online) Ratios of Fourier transforms of the 
field energy densities (integrated over k z ) at various times; 
the gauge potentials were transformed to Coulomb gauge. 



In the non-Abelian case, however, on three- 
dimensional lattices transverse magnetic fields are much 
less dominant (see, e.g. Fig. 5 in [l2|) although they do 
form larger coherent domains in the transverse plane at 
intermediate times than E± , Fig. [6l Longitudinal fields 
and locally non-zero Chern-Simons number ~ tr E • B 
emerge, also. Nevertheless, Fig. [7] shows that E z > B Zl 
aside from B± > E±. Hence, the field configurations 
are such that particles are deflected preferentially in the 
longitudinal z-direction (to restore isotropy). 

A third contribution to p z broadening in an expanding 
plasma, not considered explicitly here, is due to a lon- 
gitudinal collective flow field which "blows" the jet frag- 
ments to the side [17|. This mechanism is also available 
for collision dominated plasmas with (nearly) isotropic 
momentum distribution. However, rather strong flow 
gradients seem to be required to reproduce the observed 



FIG. 7: (Color online) Ratios of field energy densities. 



broadening of midrapidity jets (the flow velocity has to 
vary substantially within the narrow jet cone). In con- 
trast, color fields will naturally deflect particles with 
lower momentum by larger angles (Ap <~ E,B): the jet 
profile broadens even if the induced radiation is exactly 
collinear. It is therefore important to determine, experi- 
mentally, whether the asymmetric broadening is related 
to the macroscopic collective flow or to an anisotropy of 
the plasma in the local rest frame. More detailed sim- 
ulations should account, also, for the fact that small- a; 
gluons are already correlated over large rapidity inter- 
vals at the initial time I18I1. 
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